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Type la supemovae (SNe) serve as a fundamental pillar of modern 
cosmology, owing to their large luminosity and a well-defined re- 
lationship between light-curve shape and peak brightness ^ ^. The 
precision distance measurements enabled by SNe la first revealed 
the accelerating expansion of the universe^ now widely believed 
(though hardly understood) to require the presence of a myste- 
rious "dark" energy. General consensus holds that Type la SNe 
result from thermonuclear explosions of a white dwarf (WD) in a 
binary system^ however, little is known of the precise nature of 
the companion star and the physical properties of the progenitor 
system. Here we make use of extensive historical imaging obtained 
at the location of SN 2011fe/PTFllkly, the closest SN la discov- 
ered in the digital imaging era, to constrain the visible-light lumi- 
nosity of the progenitor to be 10-100 times fainter than previous 



limits on other SN la progenitors. This directly rules out lumi- 
nous red giants and the vast majority of helium stars as the mass- 
donating companion to the exploding white dwarf. Any evolved 
red companion must have been born with mass less than 3.5 times 
the mass of the Sun. These observations favour a scenario where 
the exploding WD of SN2011fe/PTFllkly accreted matter either 
from another WD, or by Roche-lobe overflow from a subgiant or 
main-sequence companion star. 

SN 201 Ife/PTFl Ikly (hereafter SN 2011fe) was discovered in the 
Pinwheel Galaxy (Messier 101 ; MlOl), a spiral galaxy at a distanceEHl 
of 6.4 Megaparsec (Mpc), by the Palomar Transient Factory (PTF) 
collaboration on 2011 Aug. 24 UX''. The galaxy was intensively 
monitored by our collaboration over the past decade and the site of 
SN 20 life was fortuitously observed prior to the explosion by the 
Hubble Space Telescope (HST) on several occasions. Together, these 
archival data offer a unique opportunity to constrain the nature of the 
progenitor system of SN 201 Ife. 

The absence of hydrogen and helium, coupled with the presence 
of silicon in the spectra of SNe la, are the distinguishing signatures of 
the class and exclude most massive stars as SN la progenitors. This, 
taken with the significant inferred amount of synthesised radioactive 
^^Ni (0.1-0.9 Mg)^^^ created in the explosions, suggests that the ma- 
jority of SNe la arise from rapid thermonuclear fusion of a carbon- 
oxygen white dwarf (WD) that is more massive than that made from 
single-star evolution (greater than 1 solar mass). A more massive WD 
also has the higher densities needed to ignite the runaway carbon fu- 
sion reaction and trigger the observed nucleosynthesis. Thus, a com- 
panion star that donates mass to the WD is required, even though fun- 
damental questions remain as to how the accretion of matter leads to 
the explosion for each progenitor model. 

The broadly diverse progenitor models are organised into t wo 
classes: double degenerate (DD) and single degenera te (S D)^^^^. A 
DD model involves two WDs in a close binary systeml^l^^, due to the 
release of gravitational radiation, the orbit shrinks. Ultimately, the 
lighter object (secondary) is disrupted and accretion onto the primary 
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ig nites runaway tiiermonuclear fusion. Alternatively, in SD mod- 
the primary WD accretes material from a stellar companion 
until its mass approaches the Chandrasekhar mass, Mch ~ 1.4 M0, 
at which point again a thermonuclear explosion ensues. Different SD 
models can be distinguished based on the nature of the secon dary : 
the WD can accrete either from a wind (the "symbiotic channel'l^^, 
by Roche-lobe overflow (the "RLOF channel'ii^, or by mass transfer 
from a helium star (the "helium star channel"^ '^). The secondary star 
in the symbiotic channel is often a red-giant star, while it is a subgiant 
or main-sequence star in the RLOF channel. 

Because WDs are faint and can only be observed directly in our 
own Milky Way and several very nearby galaxies, much of the focus 
in SN la progenitor searches has been on determining the nature of the 
companion star. Previous efforts (see Supplementary Table 1) have 
achieved an absolute magnitude in the V band of My ~ — 5 mag, a limit 
that is not deep enough to rule out any plausible progenitor systems. 

We have analysed the available HST observations of MlOl (see 
Supplementary Information for details). To pinpoint the precise SN 
location in the HST images, we obtained a mosaic image of the field 
of SN 201 Ife with the Near-Infrared Camera 2 (NIRC2) mounted be- 
hind the adaptive optics (AO) system on the Keck II telescope The 
Keck AO image was astrometrically registered to the HST/ ACS im- 
ages, yielding aid precision of 0".02I (or 21 mas) for the SN po- 
sition (see Supplementary Information). Figure I shows the site of 
SN 201 Ife on different scales. No object is detected at the nominal 
SN location in 4 different HST filters, or within the ~ 8tT error radius. 

This non-detection can be directly translated into limits on the 
progenitor system brightness and temperature. Indeed, each of the 
four HST bands place different limits on the progenitor brightness de- 
pending on the (unknown) effective temperature and luminosity class. 
These limits in turn can be converted to an equivalent (/-band mag- 
nitude, allowing us to quote the most restrictive single-frame image 
limit or a limit derived from a temperature-dependent stack of all im- 
ages; these quantities, as well as the absolute-magnitude limits at the 
distance of MlOl, are given in the Supplementary Information. 

Figure 2 shows the region in effective temperature vs. absolute 
V-band magnitude excluded by the HST imaging analysis. At 3000 K, 
progenitor systems of SN 201 Ife are excluded for My < I mag, and 
for temperatures larger than Tgff ~ 5000 K, the exclusion is My < —0.5 
mag. These limits rule out a symbiotic binary progenitor (a red giant 
and high-mass WD) such as RS Oph and probably T CrB^". Simi- 
larly, the helium-star binary V445 Pup^' is only marginally consistent 
with the upper limits, although the entire He-star channel cannot be 
completely excluded (see the blue-shaded area in Figure 2). Among 
the SD models, only the RLOF channel (such as U Scd^ can easily 
be reconciled with the HST constraints. Finally, all variations of the 
DD model are consistent with the non-detection of a source at the po- 
sition of SN 201 Ife, as it does not predict a bright source in the visible 
range. The brightness limits we have deduced also rule out a globular 
cluster, or any open star cluster with more than ~ 300 members as the 
site of SN 201 Ife. 

Figure 2 also shows that a mass donor with an effective tempera- 
ture less than 4800 K (spectral type redder than --^ G5) needs to have a 
zero-age main sequence (ZAMS) mass less than 2.2 M0. For stars 
with spectral type redder than an AO star (effective temperature ~ 
10,000 K), the companion star would need to be less massive than 
3.5 Mq. For Teff ^ 3000-4000 K, as expected for the red-giant branch 
(RGB) stars, the My limit excludes progenitors brighter than an ab- 
solute /-band mag Mj ~ —2. This limit is 2 mag fainter than the ob- 
served- tip of the RGB in MlOI and places an upper bound to the 
radius of if < 60 R© for Tgff = 3500K on any RGB progenitor. In a 
progenitor model that requires RLOF, this limit then demands an or- 



bital period smaller than Pg^^, = 260 to 130 days in a binary system 
with a 1.3 M0 WD (where the range of P^jb accommodates the 0.5- 
2.5 Mq range allowed for an RGB star). 

RS Oph, T CrB, and U Sco are recurrent novae in the Milky Way. 
The recurrence time of WD binaries where the WDs are close to the 
Chandrasekhar mass is expected?^- to be 10-20 yr. Though our his- 
torical images had sufficient sensitivity to detect classical novae, we 
find no evidence for any such outburst in the past 12yr at the site of 
SN201 Ife. However, we estimate a '-^37% chance that a typical nova 
could have occurred in the past --^5 yr and have been missed given the 
particular cadence of the imaging (see Supplementary Information). 

Historical imaging at other wavebands complements these visible- 
light progenitor system constraints. We have analysed 11 epochs of 
archival Chandra X-ray observations of MlOl taken in 2004 (see Sup- 
plementary Information), and derived upper limits for the X-ray lumi- 
nosity at the location of SN 201 Ife in the range (4-25) x lO-'^ erg s"' 
(depending on the details of the assumed spectrum). SD progeni- 
tor systems are thought to undergo a prolonged period (A/ ~ 10* yr) 
of steady nuclear burning during the mass-transfer process. Such 
systems should appear as luminous [Lx ~ lO-^^-IO-^^ erg s^' X-ray 
(kT ~ 100 eV)] sources, and indeed nearly 100 of these super-soft 
sources (SSSs) have been identified to da te in the Milky Way and 
other nearby galaxies, including MlOl itseli^^^^DD progenitor sys- 
tems have also been predicted to emit X-raysl2322l vvith Lx ~ 10^*- 
10^^ erg s^' . Our historical X-ray limits are, unfortunately, not deep 
enough to rule out either channel. 

We have also analysed pre-explosion Spitzer mid-infrared (IR) 
data of MlOI taken in 2004 (see Supplementary Information). No 
point source was detected at the SN location, to a mid-IR luminos- 
ity of < (1-13) X lO-'* ergs^'. These limits rule out bright red-giant 
mass donors of SN 201 Ife, in the symbiotic channel, but are consis- 
tent with a relatively faint secondary star in the RLOF channel. For 
a DD merger to produce a SN la, the merger product may need to 
evolve and cool slowlyEHfor > 1 x 10* yr. This requires the radiated 
luminosity (mostly in the IR) to be < 4 x lO"'* ergs^'. The Spitzer 
non-detection limits are certainly consistent with the requirement of 
low IR luminosity for the DD models. 

One test to determine whether SN 201 Ife comes from the RLOF 
channel of the SD model, or from a DD merger, is to detect the surviv- 
ing companion in the RLOF channel. But with our current resources, 
this is a very challenging observation for a subgiant mass donor, and 
nearly impossible for a main-sequence companion, unless there is ex- 
cessive heating caused by the interaction between the companion star 
and the SN ejecta. 

The observations of SN 201 Ife and historical imaging at the SN 
site have directly excluded symbiotic binary and probably helium-star 
progenitors. This, along with the X-ray and mid-IR limits, suggest that 
the WD exploded as SN 201 Ife was accreting matter either from an- 
other WD, or by RLOF from a subgiant or main-sequence star. Given 
the observed diversity of SNe la, however, the possibility of more than 
one progenitor for the class remains. As the trove of high-resolution 
wide-field imaging grows, covering more and more future SN explo- 
sion sites, similar analyses will eventually shed light on this possibil- 
ity. 
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Figure 1 | The site of SN2011fe in IVIessier 101 as imaged by HST/ACS. The left panel is a full-view colour picture of the face-on 
spiral galaxy Ml 01 (18' x 18' field of view) constructed from the three-colour HST/ACS images taken at multiple mosaic pointings (from 
http://hubbiesite.orgi. North is up and east to the left. M101 displays several well-defined spiral arms. With a diameter of 1 70 
thousand light years across, Ml 01 is nearly twice the size of our Milky Way Galaxy, and is estimated to contain at least one trillion stars. 
The middle panel is a cutout section (3' x 3') of the left panel, centred on the SN location. SN 2011 fe is spatially projected on a prominent 
spiral arm. The right panel is a section of 2" x 2" centred on the SN location, which is marked by two circles. The smaller circle has a 
radius of our Icr astrometric uncertainty (21 mas), while the bigger circle has a radius of 9 times that. No object is detected at the nominal 
SN location, or within the 8<7 error radius. Two nearby, but unrelated, red sources are labeled as "Star 1 " and "Star 2," and are displaced 
from our nominal SN location by ~ 9(t, formally excluded as viable candidate objects involved in the progenitor system of SN2011fe. 
Credit for the left panel colour picture: NASA, ESA, K. Kuntz (JHU), F. Bresolin (University of Hawaii), J. Trauger (Jet Propulsion Lab), J. 
Mould (NOAO), Y.-H. Chu (University of III inois, Urbana), and STScl. Note: This is a reduced-sizs figure for arxiv posting. 
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Figure 2 Progenitor system constraints in a Hertzsprung-Russeil (H- 
R) diagram compared to some proposed single-degenerate progeni- 
tors. The thick yellow line is the 2a limit in V-band absolute magnitude (My) 
against effective temperature at the SN location (see text) from a combina- 
tion of the four HST filters, weighted using synthetic colours of redshifted 
stellar spectra at solar metallicity for that temperature and luminosity class. 
A more conservative limit comes from taking the single filter that most con- 
strains the stellar type and luminosity class; shown is the 2a limit assuming 
the adopted distance modulus''^ of 29.04 mag (middle light yellow curve) 
with a total uncertainty of 0.23 mag (top/bottom light yellow curve). Depicted 
are the theoretical estimates (He-star channeM^) and observed candidate 
systems (V445 Pup^^, RS Oph^", U Sco^^^^, and T CrB^O). Also plotted 
are theoretical evolutionary tracks (from 1 IVIyr to 13 Gyr) of isolated stars for 
a range of masses for solar metallicity; note that the limits on the progenitor 
mass of SN 2011 fe under the supersolar metallicity assumption are similar to 
those represented here. The grey curve at top is the limit inferred from HST 
analysis of SN2006dd, representative of the other nearby SN la progenitor 
limits (see Supplementary Information). For the helium-star channel, bolo- 
metric luminosity corrections to the V band are adopted based on effective 
temperature^". The foreground Galactic and IVIIOI extinction due to dust 
is negligible^ and taken to be Ay =0 mag here. Had a source at the 2.0 <t 
photometric level been detected in the HST images at the precise location of 
the SN, we would have been able to rule out the null hypothesis of no signif- 
icant progenitor with 95% confidence. As such, we use the 2a photometric 
uncertainties in quoting the brightness limits on the progenitor system. 



6 



Constraints on the Progenitor 
of the Type la Supernova 
PTF11kly/SN2011fe. 

Supplementary Methods 

Prior limits on SN la progenitors from optical data 

Nine SNe la with preexisting HST data of their host galaxy have been close enough (within 25 Mpc) to search for a 
progenitor (Supplementary Table 1). No progenitor has been found; only upper Umits have been possible Limits 
range from Mg = —3.9 mag in the case of SN 2004W to Mj = —8.3 mag for SN 2003cg ^ . These limits provided only 
poor constraints — in the case of SN2006dd and SN 2006mr in NGC 1316, ruling out normal stars with initial masses 
greater than 6 M© at the tip of the asymptotic giant branch (AGB), young post-AGB stars with initial masses greater 
than 4 M©, and post-red- giant stars with initial masses greater than 9 M© ^. 

HST data reduction 

Search for the progenitor system of SN2011fe 

We searched the HST data archive and found six datasets that contain the site of SN2011fe. These include the ACS 
images observed in the F435W, F555W, and F814W filters in program GO-9490 (PI: K. Kuntz) on 2002 Nov. 13, a 
pair of images observed with the Wide Field Planetary Camera 2 (WFPC2) in the F606W filter in program GO-7909 
(PI: S. Casertano) on 1998 May 22, and a Wide Field Camera 3 (WFC3) image observed in the F469N filter in program 
GO-11635 (PI: M. Shara) on 2010 Apr. 7. More details on the available HST AaiSi are listed in Supplementary Table 
2. 



To initially locate the site of PTFl Ikly in these images, we utilised the PTF discovery images reported by Nugent et 
al. (201 la) ^. We followed the general procedure to identify SN progenitors in high-resolution //primages as reported 
by, for example, Li et al.(2005)^. We first registered the PTF images to a deep r-band MegaCam image downloaded 
from the 3.6-m Canada-France-Hawaii Telescope (CFHT) science data archive'. 12 stars were used in the astrometric 
solution for this step, yielding a precision of la = 0".055. The CFHT MegaCam image was further registered to the 
ACS/F814W image. The astrometric solution, using 16 stars, has a precision of 1(7 = 0".017. The 1(7 uncertainty of 
the final SN location in the ACS/F814W image is 0".058 = 58 mas. No object is visible at the nominal position of 
SN2011fe in any of the 7/5r images. This finding was reported as an astronomer's telegram (ATEL; Li et al. 201 la) ^. 
We noted in our ATEL the presence of two red sources within roughly the 2 x la uncertainty of our solution. Further 
analysis of these two sources was reported by Li et al. (201 lb) ^. 

To improve our astrometric solution, on 201 1 August 26 we observed the field of PTFl Ikly with the Near-Infrared 
Camera 2 (NIRC2) mounted behind the adaptive optics (AO) system on the Keck 11 telescope^. Beginning at 5:33, 
we obtained a series of H-band exposures, each consisting of 30 nondestructive readouts of 2 s each, for a total time 
on source of 360 s. The SN was sufficiently bright to use as the guide star for the tip-tilt correction. Images were 
reduced using standard IRAF^ routines, using a median combination of the (nonaligned) dithered images to correct 

' This research used the facilities of the Canadian Astronomy Data Centre operated by the National Research Council of Canada with the support 
of the Canadian Space Agency. 

^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in 
Astronomy (AURA), Inc., under cooperative agreement with the National Science Foundation (NSF). 
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Supernova 


Distance 
(Mpc) 






Limits (mag)" 






Reference 


SN 2003cg 


17.3 


Mr = 


-7.1 


M/ = 


-8.3 








SN 2003gs 


17.9 


Mv = 


-6.2 












SN 2003hv 


21.3 


Mv = 


-6.3 


M/ = 


-6.8 








SN 2004W 


17.3 


Mg = - 


-3.9 


Mz = 


-5.5 








SN 2005df 


24 


Mb = 


-6.2 


Mr = 


-5.1 








SN 2006dd^' 


19 


Mb* = 


-5.4 


Mv = 


= -5.3 


M] = 


-6.0 


2 


SN 2006mr<^ 


19 


Mb* = 


-5.6 


Mv = 


= -5.4 


M,= 


-6.0 


2 


SN 2007sr 


12-25 


Mv = 


-5.5 to -4.1 


Mi = 


-7.5 to -5.6 






3 


SN 2007on 


20 


Mg = - 


-4.5 


Mi = 


-4.5 
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Supplementary Table 1 — Limits on SN la progenitors from preexisting HST 6a\.a. "Filter Key: B = F450W, 
B* = F435W, g = F475W, V = 555W, R = F606W, / = F814W, Z = F850LP. Photometric system is AB for 
ref. 2, Vega for ref. 4. Others are not stated. ''SN2004ab at 14.1 Mpc has no constraints due to the high 
surface brightness in the galactic bulge. 'Both SN 2006dd and SN 2006mr were located in the same galaxy 
NGC 1316. 



Dataset 


Date(obs) 


Exptime(s) 


Instrument 


Aperture 


Filter 


Prop. No. 


Label" 


J8D6A2011 


2002-11-13 


900.0 


ACS 


WFC 


F435W 


9490 


ACS/F435W 


J8D6A2021 


2002-11-13 


720.0 


ACS 


WFC 


F555W 


9490 


ACS/F555W 


J8D6A2031 


2002-11-13 


720.0 


ACS 


WFC 


F814W 


9490 


ACS/F814W 


U4K2KF02R 


1998-05-22 


260.0 


WFPC2 


WFALL 


F606W 


7909 


WFPC2/F606W 


U4K2KF01R 


1998-05-22 


260.0 


WFPC2 


WFALL 


F606W 


7909 


WFPC2/F606W 


IB3P17010 


2010-04-08 


6106.0 


WFC3 


UVIS 


F469N 


11635 


WFC3/F469N 



Supplementary Table 2 — Available HST6ata at the location of PTF1 1 kly/SN 201 Ife. 
"These labels are used in the text to represent the images. 
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Supplementary Figure 1 — Astrometric solution using tine Keck AO image. Tine Keck AO image is sinown 
in tine rigint panel, and is astrometrically matched to a section of //5r ACS/F814W image on the left panel. 
The field of view is 1 ' x 1 '. Detected objects in the Keck AO images (and the counterparts in the HST image) 
are marked with black squares. 13 of them with reasonably accurate positions are used in the astrometric 
solution, achieving a la precision of 21 mas. The red square is centred on the SN and has a scale of 2" x 
2", and matches the scale of the right panel presented in Figure 1 of the main text. The green lines mark 
the boundary of the AO image in the ACS/F814W image. 

for the sky background. Prior to registration, we applied the distortion correction provided by the Keck Observatory^. 
The dithered images were registered and stacked to produce a single image with a resolution of 0".04 pixel^^ slightly 
superior to the resolution of the ACS images (0".05 pixel-i). 

We were able to identify 13 stars around the SN position that are present in both the Keck AO and ACS/F814W 
images (Supplementary Figure 1). The astrometric solution of this set of objects is limited by the difficulties of 
accurately measuring the centroid positions for the stars on the Keck AO images due to the signal-to-noise ra- 
tios (S/N) of the detections. Nontheless, the solution has a precision of Id = 0".021 or 21 mas, much improved 
over our reported precision (58 mas) in Li et al. (201 la)^. The SN location in the ACS/F814W image (filename: 
HST_9490_a2_ACS_WFC_F814W_sci.fits) is X = 3309.30, Y = 1750.12, with an uncertainty of 0.30 pixel in both 
axes. This position is consistent with our preliminary analysis^ to within la uncertainty, and also agrees with the 
position reported by Smartt et al. (2011)'". The precise position for the SN location is a = 14''03'"05'.733, 5 = 
H-54°16'25".18 (J2000.0) in the world coordinate system (WCS) of the ACS/F814W image. 

The site of PTFllkly in the other HST data reported in Supplementary Table 2 was carefully inspected, and no 
object was detected in any of the images. In the next section, we derive the limiting magnitudes of the HST data. 



^See http://www2.keck.hawaii.edu/inst/nirc2/forReDoc/post_observing/dewaip. 



Label 



3.0a 2.5a 2.0a 



ACS/F435W 
ACS/F555W 

ACS/F814W 



27.42 
27.04 
26.35 
26.60 



27.62 
27.25 
26.55 
26.80 



27.87 
27.49 
26.81 
27.05 



WFPC2/F606W 



Supplementary Table 3 — Limiting magnitudes of the HST images in the Vega magnitude system. 
Limiting magnitudes 

We employed several methods to derive limiting magnitudes for the HST images. For method I, the photometry 
for many stars near the SN position was measured, and the measured magnitude that gives an error of 0.2 mag is 
considered to be the 5 a upper hmit. The 3 a upper hmit is then 0.4 mag deeper. To perform photometry of the stars, 
we have followed the recipes described by Sirianni et al. 2005 for ACS and by Holtzman et al. (1995a, 1995b) ^^'^^ 
for WFPC2. For method II, we measured the variance of the sky background in the immediate vicinity of the SN. 
This variance, along with the instrument characteristics, was used to calculate the brightness for objects with different 
S/N. For method in, we made use of the Exposure Time Calculator (ETC) available from STScI. Assuming a flat 
continuum in Fx, we derived the magnitude for an object that would be detected at different S/N with the setup of the 
HST observations. The three methods yield limiting magnitudes that are in general within 0.1 mag of each other, and 
their average values are reported in Supplementary Table 3. 

While the limiting magnitudes of the individual HST datasets already reach great depth, the ensemble of all the 
HST data, when considered together, offers even greater depth. We focus this analysis on the V band. Since the HST 
data were observed with different filters, they can only be combined if the spectral type of the object of interest is 
known. We have made use of the spectrum database for the different spectral types published by Castelli & Kurucz 
(2004) augmented with blackbody spectra for extremely hot stars, and Pickles (1998) spectra of very cool stars 
and brown dwarfs. 

For each spectral type, we first derived the tightest limiting magnitude for the V band from the 4 individual 
HST images. Using the "synphot" package in IRAF/STSDAS, we converted the hmiting magnitudes derived for 
the ACS/F435W (B band) and ACS/F814W (/ band) into effective hmiting magnitudes in the V band. The deepest 
V-band magnitude among the 4 measurements is then chosen as our "1-frame" limiting magnitude in the V band. 
Not surprisingly, when the spectral type is blue (such as O and B stars), the HST data observed with the blue filter 
(ACS/F435W) provide the deepest V-band limit, while for the red spectral types (K and M stars), the HST data in 
ACS/F814W are the deepest. 

We then derived the hmiting magnitude for the stacked HST data. We employed two methods to do this. For the 
first method ("time- stacked" method, hereafter), we converted the HST data observed with ACS/F435 W, ACS/F8 14W, 
and WFPC2/F606W into "effective" ACS/F555W observations. We discuss ACS/F435W as an example. For a given 
spectral type, the B-band hmiting magnitude was converted to the V band using synphot. Using ETC, we derived 
the "effective" exposure time at which ACS/F555W would reach the same V-band limiting magnitude (for the said 
spectral type). The stacked //5rdata were then considered as an ACS/F555W observation with an exposure time that 
is the sum of all the "effective" exposure times, and the limiting magnitudes were derived using ETC. In general, the 
linniting magnitude for the stacked image is deeper than the "1-frame" hmits by 0.2-0.8 mag. 

We also used an altemate method ("flux-stacked" method, hereafter) to determine an effective limiting V magnitude 
from "stacking" the HST images. We began with the same synphot calculated colours for the full range of spectral 
types used above. We then calculated the goodness-of-fit statistic for each spectral type, varying the V magnitude 
over a wide range, and comparing the model flux with the measurements. The calculation was performed in flux 
units, with the uncertainty in the measured flux derived from observed 3 a upper limits given in Supplementary Table 
3. We then tabulated 2a, 2.5a, and 3a effective limits to the V magnitude by determining where the fit probabihty 
reaches 0.9545, 0.9876, and 0.9973, respectively. The method is illustrated in Supplementary Figure 2, which shows 
the calculation for a KO III model. Our results with this method were similar to those using the effective exposure time 
method above, and so we have averaged the two in determining our final constraints. 
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Supplementary Figure 2 — Illustration of the method to determine the effective limiting y magnitude 
from "stacking" the HST observations. The thick arrows show the observed upper limits in ACS/F435W, 
ACS/F555W, WFPC2/F606W, and ACS/F814W from left to right (all in the VEGAMAG system), and are 
hatched at 3c7, 2a, and la (in flux units) from top to bottom. The open diamonds (connected by dashed 
lines) show the synphot-calculated magnitudes for a KO III star over a range of V brightnesses. We calcu- 
lated the %^ statistic as a function of V magnitude to determine the effective y magnitude limits from all the 
data combined. In this example we find y > 28.19, 27.95, and 27.75 mag at 2a, 2.5a, and 3a, respectively 
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3.0(7 2.5(7 2.0(7 



Sp. Type 


"r'-frame 


stacked 


"1 "-frame 


stacked 


"l"-frame 


stacked 


03 V 


27.74 


28.02 


27.94 


28.22 


28.19 


28.47 


05 V 


27.74 


28.03 


27.94 


28.23 


28.19 


28.48 


BO V 


27.70 


28.00 


27.90 


28.20 


28.15 


28.45 


B5 V 


27.58 


27.93 


27.78 


28.13 


28.03 


28.38 


AOV 


27.44 


27.86 


27.64 


28.06 


27.89 


28.31 


A5 V 


27.31 


27.82 


27.51 


28.02 


27.76 


28.27 


FO V 


27.15 


27.78 


27.35 


27.98 


27.60 


28.23 


F5 V 


27.10 


27.78 


27.30 


27.98 


27.55 


28.23 


GOV 


27.10 


27.78 


27.30 


27.98 


27.55 


28.23 


G5 V 


27.10 


27.78 


27.30 


27.98 


27.55 


28.23 


KO V 


27.21 


27.81 


27.41 


28.01 


27.66 


28.26 


K5 V 


27.77 


28.05 


27.97 


28.25 


28.22 


28.50 


MOV 


28.31 


28.41 


28.51 


28.61 


28.76 


28.86 


M4 V 


— 


29.25 


— 


29.45 


— 


29.70 


1^15^ 


— 


29.61 


— 


29.81 


— 


30.06 


BO III 


27.69 


28.00 


27.89 


28.20 


28.14 


28.45 


B5ni 


27.59 


27.93 


27.79 


28.13 


28.04 


28.38 


Gom 


27.10 


27.78 


27.30 


27.98 


27.55 


28.23 


G5 III 


27.20 


27.80 


27.40 


28.00 


27.65 


28.25 


KO III 


27.39 


27.86 


27.59 


28.06 


27.84 


28.31 


K5 III 


27.89 


28.12 


28.09 


28.32 


28.34 


28.57 


MO III 


28.25 


28.36 


28.45 


28.56 


28.70 


28.81 


M5 IIP 


— 


29.28 


— 


29.48 


— 


29.73 


MIO IIP 


— 


30.89 


— 


31.09 


— 


31.34 


BO I 


27.67 


27.98 


27.87 


28.18 


28.12 


28.43 


B5I 


27.57 


27.93 


27.77 


28.13 


28.02 


28.38 


AO I 


27.48 


27.89 


27.68 


28.09 


27.93 


28.34 


A5I 


27.44 


27.87 


27.64 


28.07 


27.89 


28.32 


FOI 


27.36 


27.84 


27.56 


28.04 


27.81 


28.29 


F5I 


27.23 


27.81 


27.43 


28.01 


27.68 


28.26 


GO I 


27.11 


27.78 


27.31 


27.98 


27.56 


28.23 


G5I 


27.40 


27.85 


27.60 


28.05 


27.85 


28.30 


KOI 


27.51 


27.90 


27.71 


28.10 


27.96 


28.35 


K5I 


28.20 


28.32 


28.40 


28.52 


28.65 


28.77 


MO I 


28.19 


28.32 


28.39 


28.52 


28.64 


28.77 


M2I 


28.77 


28.78 


28.97 


28.98 


29.22 


29.23 


BBl'' 




27.91 




28.11 




28.36 


BB2'' 




27.94 




28.14 




28.39 


BB3* 




27.96 




28.16 




28.41 



Supplementary Table 4 — Limiting magnitudes in tine V band for tine "stacked" HST image, "fox these 
spectral types, the spectra from Pickles (1998)''^ were used, and only the "flux-stacked" limits were derived. 
^Biackbody spectra with Teg = 35,000 K (BB1), 60,000 K (BB2), and 100,000 K (BBS). For these spectra, 
only the "flux-stacked" limits were derived. 
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The reddening toward SN 201 Ife is negligible in the Galaxy \E{B - V) = 0.009 mag] and in MlOl (Nugent et al. 
2011b). Adopting a distance modulus of n = 29.05 ±0.23 mag'^'^*^, the V-band magnitude limits in Supplementary 
Table 3 can be converted to absolute magnitudes M{V). We fit a cubic spline function to the 2a limiting magnitudes as 
a function of the effective temperatures of the different spectral types, and present the results in a Hertzspnmg-RusseU 
diagram (H-R diagram, Figure 2 of the main text). 

For a given effective temperature and My limit the minimum effective radius /?eff,imn of the allowed progenitor can 
be tabulated assuming a "bolometric correction" from the V-band to the total bolometric luminosity [BCv(reff)] At 
Teff = 3000 K, BCv = -5.00 mag; at Jeff = 4000 K, BCy ==-1.12 mag. We take /?eff,min to be 

^q8.4694-0.2[BCv {T^b)+Mv] 

^eff,min/^© = =2 

''eff 

where My is the limit we calculate above. At Jeff = 3010,3490,4050 K we find the Umit on My = 0.59,0.24, -0.22 
mag and /?eff,min = 240, 63, 327?0 (with Rq the mass of the Sun). 

Implications on tlie progenitor system of SN 2011fe from tlie H-R diagram 

Type la supernovae (SNe la) are believed to be caused by the thermonuclear explosion of a white dwarf composed 
mainly of carbon and oxygen (i.e., a CO white dwarf). In the most popular models, this is due to the degenerate ignition 
of carbon when the mass of the CO white-dwarf approaches the Chandrasekhar mass (~ I.4M0), the maximum 
possible mass.^ 

One of the main unsolved questions is what makes a white dwarf grow to the Chandrasekhar mass. There are 
numerous progenitor models which broadly follow into two classes: double-degenerate (DD) models^^'^^ and single- 
degenerate (SD) models^'^'^^ . In the double-degenerate models, two CO white dwarfs in a close binary are driven 
together by the emission of gravitational waves until the Ughter white dwarf is disrupted dynamically, leading to the 
merger of the two. A SN la may occur when the combined mass approaches the Chandrasekhar mass. In the single- 
degenerate models, the mass of the CO white dwarfs grows by accretion of matter from a companion star. Different 
SD models can be distinguished based on the nature of the mass donor: this can be either a main-sequence star or 
subgiant ("super-soft channel" 2*), a giant star ("symbiotic channel" 2^), a helium star ("heUum-star channel" ^'^•2^), or 
perhaps even a degenerate object. 

All of these proposed models have their advantages and problems, and there are numerous observational constraints 
ruling out particular models in individual cases (see, e.g., Podsiadlowski et al. 2008 for further discussion). In 
addition, it is quite possible, perhaps even likely, that more than one progenitor channel is required to explain the 
observed diversity among SNe la. 

Because it is relatively nearby, SN PTFllkly provides a unique opportunity to constrain the various proposed 
progenitor models, at least for this SN. Figure 2 in the main text shows the upper My limits in an H-R diagram 
from deep pre-SN HST images. The figure also indicates the location of various observed systems in our Galaxy 
that have been suggested to be prototypical SN la progenitors systems in the SD channel: RS Oph and T CrB are 
symbiotic binaries with orbital periods of 455 and 228 days, respectively, representing the symbiotic channel-'"; U Sco 
is a supersoft source with an orbital period of 30 hr^'; and V445 Pup is a binary with a heUum-star donor which 
underwent a helium nova in 2000 (its orbital period has not yet been determined) In all of these systems, the white 
dwarfs are believed to have a mass close to the Chandrasekhar mass, making them excellent candidates for SN la 
progenitors.^ 

As Figure 2 shows, the //primages already constrain possible progenitor models. In particular, it strongly argues 
against a symbiotic binary progenitor, ruling out a system like RS Oph and probably T CrB. Similarly, a helium-star 
binary such as V445 Pup is only marginally consistent with the upper limits, although it cannot be completely ruled out 
(see the blue-shaded area labelled "He-star channel" in Figure 1). Among the SD models, only the supersoft channel 

^There are also alternative models in which the exploding white dwarf is below the Chandrasekhar mass, so-called sub-Chandrasekhar models. 
In these models, carbon may be ignited by a detonation wave caused by the detonation of helium in a shell surrounding the CO core (i.e., "double- 
detonation models" ^"-^^ ). 

'As a caveat, it should be noted that it is not established in any of these systems whether the white dwarf is a CO white dwarf or an ONeMg 
white dwarf. In the latter case, the white dwarf would not undergo a thermonuclear explosion when it reaches the Chandrasekhar mass, but rather 
experience collapse to a neutron star, not producing a SN la. 
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can easily be reconciled with the HST constraints. Finally, all DD models are consistent with the non-detection of a 
source at the position of PTFl Ikly, as it does not predict a bright source in the visible range. 

RS Oph, T CrB, and U Sco are recurrent novae. Indeed, the recurrence time of WD binaries where the white dwarf 
is close to the Chandrasekhar mass is expected to be 10-20yr^-'. As Supplementary Figure 3 shows, no nova has been 
detected for the last 12yr. However, this does not rule out this channel, as the recurrence time just before the SN is 
quite uncertain and recurrent novae are only very bright for a very short period that could easily have been missed. 
Indeed, as the analysis in the following section shows, there is only a 63% probability that we would catch a nova 
outburst had it occurred randomly in the 6 yr preceding the SN explosion. 

In summary, the progenitor models that are most consistent with the historical constraint are either the supersoft 
SD channel or any DD model. 

Notes on the constraints from observed SN la progenitor candidates 

Figure 2 in the main text indicates the location in the H-R diagram of observed SN la candidates. For these 
systems we plot the effective temperature based on estimates of the spectral type (if available) or theoretical modelling 
of the donor star. The absolute magnitudes are taken from observations of these systems in quiescence and include, 
in addition to the hght from the donor star, the contributions from the hot accreting white dwarf and the accretion 
disk, which may be quite substantial. The large uncertainties in the inferred absolute magnitudes are largely due to 
uncertainties in the distance estimates and the applicable extinctions. 

RS Oph and T CrB 

RS Oph and T CrB are symbiotic binaries where the white dwarf is close to the Chandrasekhar mass and ex- 
periences recurrent novae. RS Oph during quiescence has my varying between ~ 11 and ~ 11.5 mag'"* . Recent 
distance estimates range from 1.4kpc (from modelling of the last outburst) to 3.1 kpc (assuming the secondary is 
filling its Roche lobe; see Barry et al. 2008^^ for discussion). Using an estimate of the visual extinction of 2 mag 
(lijima 2008)^^ then leads to an estimated range of My from —1.2 to —3.5. The spectral type of the donor has been 
determined to be MO to M2 III (Anupama & Mikolajewska 1999)^', implying an effective temperature of 3750 to 
3900K (Straizys & KuriUene 1981)^1 

RS Oph has a quiescent my « 10 mag. With a distance estimate in the range of 0.9 to 1.3 kpc and an adopted 
extinction of ^ 0.5 mag (Belczynski & Mikolajewska 1998)^^, this leads to an absolute My of —1 to +0 mag. Its 
estimated spectral type of M3/4 III impUes an effective temperature for the donor of 3550 to 3650 K (Straizys & 
KuriUene 1981)^^ 

USco 

U Sco is supersoft recurrent nova with a very massive white dwarf. Since its distance is not well established, 
we plot typical parameters for a system of this type from the modelling by Podsiadlowski (2003)^*^.^ Adopting a 
reasonable range in properties leads to a range in My Of 2.5 to 3.9 mag and a range of Tgff of 5200 to 6000 K. 

V 445 Pup 

V 445 F*up experienced a helium nova in 2000 and is a prototypical system for a SD model with a helium-star 
donor. Following Kato et al. (2008)^^, we adopt a quiescent My of 14.5 mag. Its distance is uncertain but has been 
constrained to be between 3.5 and 6.5 kpc, and its Ay has to be estimated to be ~ 1 .6 mag (Kato et al. 2008). Taking a 
range of 1 to 2 mag for the latter leads to a range in the absolute magnitude of V445 Pup during quiescence of — 1 .5 to 
+1 mag. The effective temperature of the donor star is not well determined, but helium stellar evolution tracks (Fig. 6 
in Kato et al. 2008) suggest a temperature betwen 25,000 K and 60,000 K. 

He-star channel 

For the He-star channel we take the theoretical models of Liu et al. (2010)^**, taking the properties of the donor 
star at the time of explosion (their Fig. 6) and applying appropriate bolometric corrections from Torres (2010) 

*These parameters are typical for progenitors in the supersoft channel, regardless of whether the white dwarf in U Sco is a CO or ONeMg white 
dwarf. 
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Historical light curves 

Regular ground-based optical imaging of MlOl has been obtained with four facilities during the course of the past 
^^13 years: the 0.76-m Katzman Automatic Imaging Telescope (KAIT) on Mount Hamilton, CA, as part of the Lick 
Observatory Supernova Search (LOSS;"*'); the Palomar 60-in telescope (P60;'*^) as part of the P60 Fast Transients in 
Nearest Galaxies survey (FasTING;^^); the 3.6-m Canada-France-Hawaii Telescope on Mauna Kea; and the Palomar 
48-in as part of the ongoing Palomar Transient Factory In addition to the ground-based imaging, MlOl was 
observed multiple times by HST; see the discussion in the previous section of the supplementary material. AH of the 
ground-based images were reduced using standard techniques, including removal of the bias and dark current from 
each frame followed by proper flatfielding. 

The projected position of PTFl Ikly/SN 201 Ife lies very close to a spiral arm structure within MlOl, leading to a 
complex background region in the vicinity of the SN (see Figure 1). Thus, in order to search for any variable optical 
emission that may be associated with the progenitor system of PTFl Ikly, we employed image subtraction to remove 
the complex, but presumably static, background signal. For each individual image (hereafter, new image), a reference 
image was constructed using a deep coaddition of several images taken with the same camera and in the same filter but 
during a different year than the new image to ensiu'e that any possible variable emission in the new frames is not also 
present in the reference image. Following the construction of appropriate reference images, for each new image the 
corresponding reference was then astrometrically aligned with the new image, both images were convolved to have 
the same seeing, and the flux in the reference image was then subtracted from the flux in the new image. 

Following proper subtraction, we examined whether any variable emission can be detected at the location of the 
SN. We measured the flux in an aperture centred on the SN position, which is precisely measured to an accuracy 
< 1 pixel by transforming a P48 image that contains the SN to the same coordinates as the subtracted image. To 
determine the uncertainty in this single-aperture measurement, we measured the variance in 1000 apertures randomly 
placed within a ^ 90" x 90" box centred on the SN position, which ensures similar background properties as the 
SN position but excludes strong systematic variations related to residuals associated with poorly subtracted bright 
stars. Using this method there are multiple images where we detect flux at the SN position with greater than 3(7 
confidence; however, since we have ^ 3000 images of MlOl, we would expect about five 3(7 detections based on 
Gaussian fluctuations. In order to test the null hypothesis that no emission is present in the subtraction images, with 
the alternative hypothesis that emission is present in an image, we employ the false-discovery rate method (FDR; 
see^^,^^) which controls for the fact that in multiple hypothesis testing, in this case many images, there are multiple 
comparisons to the null hypothesis. Using the FDR method, and setting a = 0.05 (equivalent to a 2c7 detection), we 
are unable to reject the null hypothesis that no variable emission is present at the position of the SN. This result is 
not surprising given that the typical individual images from KAIT (OTciear(3c7) ss 18.8 mag), P48 (»iMouid-R(3c7) w 21 
mag), and P60 (m^(3c7) 21 .5 mag) are too shallow to detect all but the very brightest novae (see Figure 3), and even 
that would require an image to be taken within a few days of the peak brightness of the nova. While the CFHT images 
are deep enough to detect even the very faintest novae, since only 7 images were taken over slightly more than 7 yr the 
probabiUty of a nova being bright enough for detection concurrently with an acquired image is only ~ 7%. 

To further constrain the possibility of pre-SN variability, we stack individual images grouped closely together in 
time in order to obtain limits that are deeper than those in the individual frames. Should variable emission, such as that 
associated with a nova, be present within the data, then the stacked images will reflect the average observed flux from 
each of the individual images. On week-to-month timescales the detection threshold for stacked images decreases 
faster than the average flux from a nova, and therefore stacked images will allow an improved opportunity for the 
detection of variabihty on timescales similar to those of a nova. For KAIT we stacked images on intervals of 30 days 
and for entire observing seasons. These images extend over the full lifetime of LOSS, ~ 13 yr, and are shown in black 
in Figure 3. For the P60 we employ both 7 day and 30 day stacks which extend from late 2006 through mid 2009, 
shown in green in Figure 3. For the P48 observations of MlOl we employ 7 day and 30 day stacks extending over the 
full lifetime of PTF from mid-2009 to the present. These limits are shown in red in Figure 3. Following the production 
of stacked images, we use the same procedure described above to test for variable emission: we subtract the flux from 
a template image and measure the fliK at the SN position compared with several randomly placed apertures in the 
vicinity of the SN position. Using the FDR method, we once again find no significant detection of variable emission 
at the position of the SN. 

Following the non-detection of any pre-SN variabihty, we attempt to constrain the probability that a nova could 
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have occurred in the 13 yr prior to PTFl Ikly/SN 201 Ife but that gaps in our observations or insufficient depth in our 
stacked images prevented us from detecting the signal. The true characterisation of this probabiUty requires an a priori 
knowledge of the nova luminosity function (NLF) including the distribution of nova rise times, which are typically 
short (see, e.g.,^^). To the best of our knowledge, a complete sample of the distribution of the NLF has not been 
assembled to date, and thus we must use proxies in our estimates of the probability of detecting a nova signal. As a 
proxy for the true NLF we use the light curves of two novae discovered during the FasTING survey in M31, nova 2008- 
09a and nova 2008-09c^^. We elect to use extragalactic novae to avoid the large systematic uncertainties associated 
with measuring the distance and reddening toward Galactic novae, and in particular, the two selective novae have 
strong constraints on the rise time as a result of the high-cadence observations from FasTING. We can characterise 
these two novae using the maximum-magnitude vs. rate-of-decline (MMRD) relation (for a theoretical explanation 
see, e.g.,"^^, although see also"*^), from which we learn that nova 2008-09c is a typical nova with Mg = —7.8 mag and 
a time to decline from peak by 1 mag, h = 9.1 days'*-'. Nova 2008-09a, on the other hand, is both faint, Mg = —6.8 
mag, and fast evolving, t\ = 6.3 days^^, meaning it is among the most challenging type of nova to detect. Using a 
smooth spline fit to the data for these two novae, we deredden them and shift them to the distance of M 101 to serve as 
a proxy to the true NLF. The light curves for each of these novae are shown in Figure 3: 2008-09c is in green, while 
2008-09a is in blue. 

The typical detection limits from KAIT are too shallow to detect all but the most luminous novae, and even in these 
cases the signal would often be degraded given the large number of KATT images needed to reach a detection threshold 
of m > 21.5 mag. Therefore, we focus our efforts on constraining the presence of a nova during the time of the P60 
and PTF observations, which constrains the following Umits to the past five observing seasons starting from late 2006 
until the time of the SN (see the bottom panel of Figure 3). To determine the detection probabilities, at random times 
we insert a fake nova signal into the data and determine the average flux from the nova based on the actual observation 
times of the stacked images that overlap the nova in time. If the average flux is greater than the measured detection 
threshold from the stacked image, then the nova is considered detected. We have repeated this procedure 1 million 
times for each of our example novae, and we find a detection probability of 2008-09a-like and 2008-09c-like novae of 
0.440 and 0.628, respectively. 

It is important to note, however, that these detection limits provide a lower limit to the true detection probability 
since we have adopted a stacking scheme that does not include overlapping windows (the precise scheme we use can 
be visualised using the horizontal error bars of the limits shown in Figure 3). Under this method a nova that peaks 
directly in between two stacking windows could escape detection, as each window overlaps with only a fraction of the 
nova. To prevent this from happening, an additional window, preferably centred near the time of peak brightness from 
the nova, should be created. Given the broad complexities associated with the image sequences, for example variable 
cadence and detection thresholds, the ideal method to test for the presence of a nova would be to stack all the images 
in every possible grouping on timescales shorter than those typical for a nova, ^ 30 days. Such an effort would be 
impractical for the ~ 3000 pre-SN images of MlOl discussed here, and is beyond the scope of the current paper. 

X-ray upper limits 

We analysed 11 epochs of Chandra observations of the field containing PTFllkly, all taken in 2004. Additional 
observations in which the SN location is near one of the CCD edges were removed. We detected no bright source in 
any epochs, and we additionally sub-selected the epochs capable of producing the most stringent flux limit. These 
were the 4 epochs (ObsIDs 6114, 6115, 4735, and 6118) where the SN was within 3.5' of the observation aimpoint, 
placing the source on the most sensitive ACIS-S3 detector. The total exposure in these 4 epochs is 142 ks. 

Count rates and spectra were extracted from a 4-pixel (2") radius aperture, containing 90% of the encircled photon 
flux, around the SN position for each epoch using standard tools in CIAO 4.3^. The sky background was estimated 
using an annular extraction region with an inner radius of 20 pixels and an outer radius of 50 pixels about the source 
position. Given 9 counts (0.3-8 keV) detected in the source region over an expected background of 1 1 .2 counts, we 
determine a 3a upper limit on the source count rate of 9.1 x 10~^ count 

We fit the summed spectra in ISIS* by minimizing the Cash (1979) statistic"*^. Assuming an absorbed power-law 

' http ://cxc.harvard.edu/ciao/index.html. 
*'http://space.mit.edu/CXC/ISIS. 
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Supplementary Figure 3 — Historical limits on pre-SN emission from the progenitor system of PTFllkly/SN 
2011fe. Shown in the top panel are the limiting magnitudes for stacked images from pre-SN imaging of M1 01 
from KAIT {R band), CFHT {R/r band), P60 (g band), and RTF {Rmoum) extending > 12yr prior to the SN 
explosion. The KAIT data are stacked on both monthly (light black) and yearly (dark black) intervals, CHFT 
limits are shown in magenta, P60 data are stacked on both weekly (light green) and monthly (dark green) 
intervals, and the PTF/P48 data are stacked on weekly (light red) and monthly (dark red) intervals. The 
width of the error bars indicates the time elapsed between the first image and the last image included in a 
stack. The bottom panel shows a zoom in on the 5yr prior to the SN, when the P60 and P48 data provide 
the strongest constraints on a possible nova. Overplotted on the data are the light curves of two novae from 
M31 : 2008-09a, a quickly-fading, faint nova (blue), and 2008-09c, a typical classical nova (green;'*-'), with 
their apparent brightness scaled to how they would appear in M101 . The limits show we would have been 
capable of detecting such outbursts, with the probability of detection 44% and 63% for 2008-09a-like 
novae and 2008-09c-like novae, respectively. 
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Supplementary Figure 4 — The limit on the progenitor radius as a function of temperature assuming a 
blackbody model, derived from the X-ray luminosity limits. 

model with photon index F = 2 and Galactic column density A^h = 18 x 10^*^ cm"^^*^, we determine a 3(7 limiting 
unabsorbed flux of 4.1 x 10^^ erg s ' (0.3-8 keV). For a 67 eV blackbody (e.g., Gilfanov & Bogdan, arXiv 1103.3659; 
see, however, Izumi et al. 2010^', the bolometric flux limit is 2.5 x 10^^ erg s^^ Supplementary Figure 4 displays the 
limit on the progenitor radius as a function of temperature, assuming a blackbody model. 

Spitzer data reduction 

We analysed pre-explosion Spitzer/IRAC data taken on UT 2004 March 8.09 and 2004 March 8.49 (PID:60, PI: G. 
Reike). Using 12 stars on the FTP discovery images^ that are also present in the Spitzer images, we were able to locate 
the SN position with a la uncertainty of 0".132 (132 mas). No point source was detected in any of the Spitzer data. 
Limiting magnitudes were derived in a 20 x 20 pixel box centred on the location of PTFl Ikly; see Supplementary 
Table 4. 

According to a mid-infrared (IR) H-R diagram presented by Ardila et al. (2010)^^, the mid-lR luminosities for red 
giants (KO HI to M6 HI) are in the range of 4 x 10^^ to 4 x 10^^ erg s~^ Thus, the Umits from Spitzer data alone rule 
out the brightest red giants, but not red giants with moderate to low luminosities. 

The Spitzer mid-IR upper limits are potentially useful as constraints on the DD models. Shen et al. (2011)^^ 
suggests that, for a DD merger to produce a SN la, it needs to evolve and cool slowly, for > 1 x 10^ yr. This means 
the merger product would radiate at smaller than a fraction (^ 1 %) of the Eddington limit (^ 4 x 10^** erg s^'), i.e., < 
~ 4 X lO^^ergs"^ The nondetections in the Spitzer data are thus consistent with the requirement of low IR luminosity 
from the DD models. 
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Supplementary Table 5 — Spitzer 3a upper limits. 
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